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Abstract. In this paper we study the near-infrared emission spectrum of IRAS 16594—4656, a bipolar post-AGB 
star with spectral type B7 and no observed ionization. Using optical and near-infrared photometry we determined 
the total extinction towards this object to be Ay = 7.5 ±0.4 mag and derived a distance of 2.2 ±0.4 kpc, assuming 
a luminosity of 10 4 Lq. The near-infrared spectrum shows strong H2 emission lines and some typical metastable 
shock excited lines such as [Fell] 1.257 & 1.644 /im. We determined the rotational and vibrational excitation 
temperatures, as well as the ortho-to-para ratio of the molecular hydrogen. Based on these we argue that the 
H2 emission is mainly collisionally excited. Line ratios indicate that the H2 emission originates in a ~25 kms - 
C-type shock. On the other hand, the metastable lines, and especially the [Fell] emission lines, indicate the 
presence of a ~75 kms -1 J-type shock. Hence we postulate that the H2 emission originates where the stellar wind 
(with an observed terminal velocity of ~126 kms -1 ) is tunneled through an equatorial density enhancement, 
impinging almost tangentially upon the circumstellar material. The [Fe 11] emission either occurs along the walls 
of the bipolar lobes where the transverse shock velocity would be higher, or could originate much closer to the 
central star in shocks in the post-AGB wind itself, or possibly even an accretion disk. Further high resolution 
near-infrared spectra are currently being obtained to confirm the proposed geometry and kinematics. 
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1. Introduction 

Post-AGB stars represent an important transition phase 
in the evolution of low and intermediate-mass stars, be- 
tween the asymptotic giant branch (AGB) and the plan- 
etary nebula (PN) phases. During this period, the de- 
tached circumstellar envelope of gas and dust is expanding 
away from the star. Meanwhile the star itself is increasing 
in temperature at about constant luminosity. This phase 
lasts a few thousand years depending upon the star's core 
mass. When the temperature is high enough and the star 
photo-ionizes the nebula, it has entered the PN phase. 

In spite of extensive study, the evolution from the 
AGB toward the PN stage is still poorly understood. The 
drastic changes observed in circumstellar structure and 
kinematics are particularly puzzling. During late AGB or 
early post-AGB evolutionary stages, the geometry of the 
circumstellar material changes from more or less spheri- 
cally symmetric to axially symmetric, with the result that 
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most PNe exhibit axisymmetric structures, ranging from 
elliptical to bipolar. Bipolar PNe are very likely to pos- 
sess molecular envelopes that are readily detectable in 
the near-infrared ro- vibrational lines of H2 (Kastner et al. 
1996). The available data suggest that the onset of near- 
infrared H2 emission in PNe can be traced back to the pre- 
planetary nebula (PPN) phase but not back to the AGB 
phase of evolution (Weintraub et al. 1998). These obser- 
vations suggest that further studies of H2 emission from 
PPNe may offer insight into the transition from AGB star 
to PN and from spherical to axisymmetric mass loss. 

The study of transition objects showing H2 emission 
at an early stage is crucial for understanding the hy- 
drodynamic processes shaping the nebulae. The H2 lines 
can reveal details about the physical conditions in the 
shocks associated with these hydrodynamic processes and 
thereby help constrain models of the interaction of the 
central star with the AGB remnant. Only 13 post-AGB 
stars with H2 emission have been detected in molecu- 
lar hydrogen (Garcia-Hernandez et al. 12002(1 . RAFGL 
2688 (the Egg Nebula) and RAFGL 618 are perhaps the 
best-studied examples of such transition objects and the 
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only two whose H2 emission has been extensively stud- 
ied. Both objects display molecule-rich circumstellar en- 
velopes, carbon-rich circumstellar chemistries, and dusty, 
bipolar reflection nebulosities. Whereas RAFGL 2688 har- 
bors an F-type central star and is not known to contain 
an H 11 region, RAFGL 618 possesses a B-type central star 
and does contain a compact Hn region, suggesting that 
the latter nebula has evolved further toward the PN stage. 
In this paper we study in detail the excitation mechanism 
of the H 2 emission in the PPN IRAS 16594-4656. 

IRAS 16594-4656 is classified as a post-AGB star for 
several reasons. In the IRAS color-color diagram it has 
the colors of a planetary nebula (PN) (Van de Steene & 
Pottasch ll993[) . but it has never been detected in the ra- 
dio (Van de Steene et al. l2(JL)(Ja|) . It has a large infrared 
excess due to dust with a color temperature of 173 K. 
It displays a double-peaked spectral energy distribution, 
with the peak in the mid-infrared much brighter than the 
peak in the near-infrared (Van de Steene et al. 1200 0a). 
It possesses a CO envelope with an expansion velocity 
of at least 16 kms -1 (Loup et al. 1990). The chemistry 
of IRAS 16594—4656 appears to be carbon-rich. This is 
based on the detection of unidentified IR emission features 
at 3.3, 6.2, 7.7, 8.6, 11.3, 12.6, and 13.4 /mi, as well as the 
21 /jm feature (Garci'a-Lario et al. I1999J) . all commonly 
associated with a carbon-rich chemistry. However, tenta- 
tive detections of crystalline silicate features (mostly py- 
roxenes) were also reported by Garcfa-Lario et al. J1999). 
The latter are associated with an oxygen-rich chemistry, 
which could point to the fact that IRAS 16594—4656 only 
turned carbon-rich during one of the most recent thermal 
pulses. 

The optical spectrum of IRAS 16594—4656 shows a 
spectral type B7 with significant reddening. The Balmcr 
lines of hydrogen are in emission, together with a weak 
[0 1] emission at 6300 A. The Ha emission has a P-Cygni 
type profile indicative of a stellar wind with a terminal 
velocity of approximately 126 kms -1 (Van de Steene et 
al. l2000bl) . 

In the optical the nebulosity is dominated by scat- 
tered light, not emission lines as in the case of plane- 
tary nebulae (Hrivnak et al. I1999J) . The HST image of 
IRAS 16594—4656 shows a bipolar morphology (Hrivnak 
et al. I1999[) . with the southwestern lobe tilted towards us 
at an intermediate orientation (Su et al. |2~00 It , as well 
as concentric arcs (Hrivnak et al. I2001|) . Su et al. ([2000 ) 
noted that in addition to the centro-symmetric polari- 
metric patterns, point-symmetric patterns are also seen 
in IRAS 16594—4656. Such patterns provide strong, inde- 
pendent evidence for the presence of a circumstellar disk 
or torus (hereafter referred to as an equatorial density en- 
hancement or EDE). 

Spectroscopic studies of H2 are particularly useful be- 
cause the intensity ratios of H2 lines arising from different 
excited levels provide diagnostics of the excitation mech- 
anism and physical conditions. There are two likely exci- 
tation mechanisms for H2 in post-AGB stars: UV pump- 
ing ("fluorescence") by stellar photons, and thermal (i.e., 



Fig. 1. The HST mage of IRAS 16594-4656 through fil- 
ter F814W (Hrivnak et al. ITM^ with the slit position 
indicated by the grey lines. North is to the top and East 
to the left. 



collisional) excitation by the gas. It is possible that both 
mechanisms play a role within a single post-AGB star. All 
spectroscopic (and imaging) studies of post-AGB stars to 
date have concentrated only on the strong H 2 (1,0) S(l) 
2.121 fini and H 2 (2,1) S(l) 2.247 spectral lines. Their 
ratio can in principle distinguish between the two exci- 
tation mechanisms. However, its utility as a diagnostic 
breaks down at high densities, when collisions modify the 
populations of the lowest vibrational levels (v = 1,2), 
driving them to their thermal values, even when the gas 
is fluorescently excited. Emission lines from higher vibra- 
tional levels are less affected, and can still distinguish gas 
in which UV photons are responsible for the H2 excita- 
tion. Furthermore time-dependent effects might enhance 
the intensity of UV-pumped emission soon after the on- 
set of far-UV irradiation, enabling fluorescent H2 emission 
to rival the high surface brightness typical of shocked gas. 
This makes it particularly important to study more than a 
handful of post-AGB stars, over a large spectral range, in 
order to achieve a broad, balanced picture of their molec- 
ular components. 

For the first time we present a full J-, H-, if-band spec- 
trum of a post-AGB star. The observations are presented 
in Sect. The extinction and distance are derived in 
Sect. Eland the infrared spectrum is discussed in Sect.0] 
Many H2 emission lines were detected and we consider 
them all to discuss the excitation mechanism in Sect. |SJ 
The results are discussed in Sect. Eland final conclusions 
are given in Sect.0 

2. Observations 

The near-infrared spectrum of IRAS 16594—4656 was ob- 
tained with SOFI on the NTT telescope at La Silla (ESO) 
in June 1998. The slit width was 1" and the slit orien- 
tation almost east-west. An image with the slit position 
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Fig. 2. The H 2 spectrum of IRAS 16594-4656. 

indicated is presented in Fig. ^ The SOFI pixel scale is 
0'.'29 per pixel. Between 0.95 and 1.65 /im the spectrum 
was obtained with the GB grism and has a spectral res- 
olution of 7 A per pixel. Between 1.53 and 2.52 /im the 
spectrum was obtained with the GR grism and has a res- 
olution of 10.2 A per pixel. The seeing was 0'.'7. The red 
spectrum is an average of 6 cycles of 10 DITs of 1.5 sec 
and the blue spectrum an average of 6 cycles of 10 DITs 
of 2 sec. The airmass was 1.45. 



Fig. 3. The other emission lines of IRAS 16594-4656. 

The spectra were reduced with an adapted version of 
the ESO package eclipse and further processed and an- 
alyzed in iraf. A selection of 3 standard stars, observed 
throughout the night, was used to construct the response 
curve. To obtain absolute flux tables for the standard 
stars, we constructed theoretical spectra for the standard 
stars based on NextGen atmosphere models produced 
with the stellar atmosphere code phoenix (Hauschildt et 
al. H999|l . which were subsequently convolved to the SOFI 
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Table 1. The observed and intrinsic magnitudes of 
IRAS 16594—4656. The intrinsic magnitudes are arbitrar- 
ily scaled to V — mag. 



B a 


ya 


Rc 
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It 


A* 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


m obs 16.31 


14.60 


13.27 


12.00 


9.73 


8.85 


8.20 


m intr -0.09 


0.00 


0.01 


0.06 


0.16 


0.18 


0.23 



a Hrivnak et al. ( 1999 



6 Van de Steene et al. ©000a) 

resolution and wavelength grid. Next we included an ex- 
tinction component and scaled the spectra to match the 
observed J-, H-, and JT-band flux values. These spectra 
were compared with the observed ones to derive the av- 
erage response curve. The object spectra were calibrated 
with this smoothed response curve. Comparison of the ob- 
served, calibrated spectra of the standard stars and the 
theoretical spectra provided us with the spectrum of the 
atmosphere, which was subsequently used to correct for 
telluric absorption. 

3. Extinction and distance 

IRAS 16594—4656 lies very close to the Galactic plane 
(b = — 3?3). As such, it is expected to be significantly 
reddened. Because the reddening must be accounted for in 
the analysis of the emission lines and the continuum, we 
begin this section with a determination of the reddening. 
We used BVRqIc photometry taken from Hrivnak et al. 
( 1999 ) and JHK photometry taken from Van de Steene et 
al. H2()()()ajl . We excluded the observed Z^band magnitude 
from our analysis since it is apparently contaminated by 
emission from either atomic lines or hot dust. The spectral 
type of the central star is B7 (Van de Steene et al. [2000b) 
and we assumed that the star has a low surface gravity. 
The intrinsic colors of the central star were derived from 
an ATLAS9 stellar atmosphere model (Kurucz 1994) with 
T off = 12000 K and log[c//(cms- 2 )] = 2.00. The data are 
summarized in Tabled 

We determined the total extinction (i.e., interstellar 
and circumstellar extinction combined) for each combi- 
nation of observed photometric bands (excluding V— Rc) 
such that the dereddened color would match the intrin- 
sic color. We then averaged all the measurements of Ay 
and determined the standard deviation. We used the ex- 
tinction law from Cardelli et al. (1989) and determined 
Rv by minimizing the standard deviation. The resulting 
values are: Ay = 7.5 ± 0.4 mag, with Ry = 4.2. Note 
that the uncertainty quoted here is the uncertainty in the 
average value, not the spread of the individual measure- 
ments. The galactic extinction estimate at the position of 
IRAS 16594—4656 is Ay = 6.6 mag according to Schlegel 
et al. I|1998ll . who assume an Ry = 3.1 extinction curve. 
This value is in good agreement with our determination. 
The extinction value is considerably higher than the esti- 
mate Ay w 3.2 mag given by Hrivnak et al. I|1999f) . The 
latter estimate pertains to the interstellar extinction com- 



ponent alone at an assumed distance of more than 1.3 kpc, 
and is based on the galactic extinction study by Neckel & 
Klare ljT35U)l . 

Visual inspection of the Digitized Sky Survey image 
taken in the direction of IRAS 16594—4656 reveals that 
the extinction is very patchy in the neighborhood of this 
object. Hence the use of estimates from galactic extinction 
studies should be considered very uncertain. On the other 
hand, some of the difference between the two extinction 
values could be due to significant circumstellar extinction, 
which would be supported by the large infrared excess of 
this object. However it is unclear how much of the total 
extinction is of circumstellar origin. 

With this value for the extinction and the flux cali- 
bration from the Kurucz model, wc can determine that 
the distance to IRAS 16594-4656 is (2.2 ± 0.4) l}J 2 kpc, 
where L\ is the luminosity of IRAS 16594—4656 in units of 
10 4 L Q . This value is in good agreement with the distance 
of 2.5 L X J 2 kpc derived by Su et al. H2()()l|l . based on the 
dereddened bolometric flux of IRAS 16594-4656. The lat- 
ter value may have been somewhat overestimated due to 
stellar radiation escaping through the lobes. The distance 
is also somewhat sensitive to uncertainties in the interstel- 
lar extinction estimate, although this effect is mild since 
most of the bolometric flux is emitted in the mid-infrared. 

4. Infrared lines 

4.1. Atomic hydrogen 

The Br7 line is in emission, but doesn't show a P-Cygni 
profile as we observed in the Ha line (Van de Steene et al. 
2000b). The Bi'7 flux detected in the SOFI spectrum is 
approximately 40% lower than measured from the higher 
resolution IRSPEC spectra obtained in 1995 by Van de 
Steene et al. (2000a) and in 1994 by Garcfa-Hernandez et 
al. l{2"0T)2"|l . For the IRSPEC observations the slit had a 
width of 4'.'5 and was oriented N-S and E-W, respectively, 
while in our SOFI observations the slit had a width of only 
1" and was oriented E-W. Based on the continuum flux, 
both calibrations agree. This could indicate that we missed 
some line emission due to the narrow slit width in the 
SOFI observations. However, this is not a very plausible 
explanation. The Br/y emission is formed in the post-AGB 
wind, very close to the central star. The emission region 
is much smaller than 1" in diameter, so if we missed B1-7 
emission, we would also have missed continuum emission 
from the central star. Reflection also seems a very implau- 
sible explanation, since it is not expected that dust scat- 
tering in the infrared would be efficient enough to cause a 
40% discrepancy in the observed fluxes. This leaves as the 
only alternative explanation that the strength of the Br7 
emission became less between 1994/5 and 1998. We have 
seen spectral variability in other, similar objects in our 
sample, so this possibility cannot be ruled out. Further 
observations will be needed to settle this point. 

The Pa/3 line has been observed for the first time in 
this spectrum, and is shown in Fig 13 
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Table 2. The emission lines detected in the spectrum of IRAS 16594—4656. The observed flux values F b s in the 5 
column are not corrected for extinction, but the F cor values in the last column are corrected for extinction. 
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3.17( 
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-Lay 


(1,0) SCO) 
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-13) 

-Lay 
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(1,0) Q(2) 
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-3) 
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-13) 


(1,0) Q(3) 
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-14) 


3.04( 
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1.22(-14) 
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-3): 
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1.2821 
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-12) 
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-14) 


Br7 
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a Entries such as ( — 14) stand for xlO 14 . 

b This line is a blend of the [C 1] 0.9824 & 0.9850 lines. 



4.2. Molecular hydrogen 

The observed H2 lines are shown in Fig. [2] The strongest 
H2 emission is seen from (1,0) Q(l), (1,0) S(l), and 
(1,0) Q(3). The weakest line is (2,1) S(l) and no other lines 
with v — 2 or higher are detected. 

Compared with the IRSPEC value given by Garcia- 
Hernandez et al. (20Q2J), we observe ^50% less flux for 
the H 2 (1,0)S(1) line. The flux of the faint H 2 (2,1)S(1) 
line however agrees for both observations (our detection 
has a large error margin though). This discrepancy could 
either be due to flux missed outside the SOFI slit, or to 
variability in the post-AGB wind. See also the discussion 
in the previous subsection and Sect. El 

4.3. Carbon 

The [C 1] 9824 & 9850 A lines are metastable lines often 
observed in shocks (Hollenbach & McKee 1989). Neutral 
carbon is readily ionized, so the existence of [C 1] requires 
a region sheltered from ionizing radiation, or indicates the 
lack of it. 

Both [C 1] 9824 A and 9850 A come from the same up- 
per level and have a fixed intensity ratio of approximately 
1:3. The fact that [Ci] 9824 A is not clearly detected in 
our spectrum is most likely due to telluric absorption. 

4.4. Iron 

[Fen] lines are known to be good tracers of astrophysi- 
cal shocks in supernovae, Herbig-Haro objects, PNe, etc. 
[Fc 11] lines are very sensitive to the high densities and tem- 
peratures in the shocked gas. The excitation temperature 
of the upper term of the IR multiplets of [Fe 11] , which fall 



in atmospheric windows at 1.2 and 1.6 /im, is 2-3 times 
lower than of optical lines and the critical densities of these 
[Fen] lines are ^10 times higher. 

The [Fen] a 4 F-a 4 D 1.644 /im over Br7 intensity ratio 
is often used as an indicator of shock excitation. The ratio 
expected for shock-excited gas is much larger than 1, but 
the [Fen] 1.644 fxm/Brj ratio expected for radiatively ex- 
cited gas is only approximately 0.06 (Graham et al. H987fl . 
The ratio which we obtain for the dereddened line ratios 
is 0.9, which is much larger than the value typically found 
in H II regions. 

5. Excitation of molecular hydrogen 

The excitation of molecular hydrogen has been the subject 
of several theoretical investigations which are reviewed 
by Sternberg : 1991 )t and Burton (1992). Two principal 
sources of H 2 emission are (1) photodissociation regions 
(PDRs) , where the molecules are vibrationally excited by 
far-ultraviolet (FUV) pumping or collisionally excited in 
gas heated by FUV radiation; and (2) shocked regions, in 
which the hydrogen molecules are collisionally excited in 
hot gas behind the shock front. In post-AGB stars both 
forms of excitation may occur. Observations of the rich 
H 2 ro- vibrational spectrum yield line ratios that are valu- 
able probes of the physical conditions within the emitting 
source. These allow us to determine the excitation mech- 
anism. 

Garcfa-Hernandez et al. H2002|) concluded that fluo- 
rescent H 2 emission becomes active when the central star 
reaches a temperature that corresponds to spectral type 
A. Shock excited H 2 emission was only detected in objects 
with a marked bipolar morphology, sometimes at an early 
stage during the post-AGB phase. Sometimes this emis- 



6 



Van de Steene & van Hoof: Shock emission in IRAS 16594—4656 



sion was localized in the waist and/or in specific regions of 
the bipolar lobes, where a physical interaction exists be- 
tween fast- and slow-moving material. In this section we 
investigate the excitation mechanism of H 2 in IRAS 16594- 
4656. 



5.1. Rotational and vibrational temperatures of 



Because H2 is formed from two identical particles, only 
quadrupole transitions are allowed, and transition proba- 
bilities are very low. As a consequence, the radiative life- 
times of the levels are long and all levels will be popu- 
lated thermally within a given vibrational state. If the 
relative strengths for a few of the transitions can be mea- 
sured, then this information can be used to derive the ro- 
tation temperature for that vibrational state. If collisions 
alone excite molecular hydrogen, then populations in dif- 
ferent vibrational states will also be in thermal equilibrium 
and all levels can be described by the same temperature. 
However, if the excitation mechanism includes a fluores- 
cent component, higher vibrational states will be overpop- 
ulated, thus indicating a vibrational temperature higher 
than the rotational value. Hence, thermal excitation will 
give similar rotational and vibrational temperatures, while 
fluorescence is characterized by a high vibrational temper- 
ature and a low rotational temperature (Rudy et al. 12001)1 . 

Fig. 0] provides a graphical determination of the tem- 
perature and excitation mechanism of the molecular hy- 
drogen in IRAS 16594—4656 in the standard manner (see, 
e.g., Martini et al. 11997)1 . The y-axis shows the log of the 
relative column density of molecular hydrogen in a given 
state, N(v,J), divided by the statistical weight, g(J). 
These were derived by converting the reddening corrected 
line fluxes F cor into column densities for the upper levels 
using the formula: 

Here all symbols have their usual meaning and Q is the 
solid angle covered by the observations. The transition 
probabilities A^i were taken from Martin et al. ©96). 
Since is unknown for IRAS 16594—4656, we present 
the data as ratios to N(l, 3)/g(3). The x-axis gives the 
temperature of the upper level of a given transition. 

Comparison of the population of upper levels with dif- 
ferent rotational quantum numbers (J- values), but identi- 
cal vibrational quantum numbers (v- values), provides an 
estimate for the rotational temperature. The temperature 
can be found directly from the slope of a straight line fit- 
ted to the data points. We have made a least-squares fit 
using the formula 

In y — In yg = a + bx, 

with x = E(v,J)/k, y — N(v,J)/g(J) and yo = 
N(l,3)/ g(3). The rotational temperature is simply given 
by T r = -1/6. The results are: a = 4.90 ± 0.30, b = 
(—6.93 ± 0.41) xl0~ 4 . Hence, the rotational temperature 
T r = 1440 ± 80 K. 



Fig. 4. The column densities for the various observed ro- 
vibrational levels of H2. The best-fit line based on a rota- 
tional temperature of T r = 1440 K is shown as well. 



The vibrational temperature is measured from the 
slope of a line passing through data points with differ- 
ent vibrational quantum numbers but the same rotational 
quantum number in the \n[N(v, J)/g(J)], T-plane. We de- 
termined T vib = 1820 ± 240 K. 

The fact that the two values differ only by 1.6 u 
is consistent with the assumption that fluorescence 
doesn't contribute strongly to the excitation and that in 
IRAS 16594—4656 the H 2 excitation is mainly collisional. 

5.2. The ortho-to-para ratio 

The ortho-to-para ratio of molecular hydrogen is the ratio 
of the total column density of ortho-H2 (all odd J states) 
to para-H 2 (all even J states). 

If the ro-vibrational states are predominantly popu- 
lated by collisional processes, the ortho-to-para ratio will 
be very close to 3. On the other hand, if FUV fluorescence 
plays an important role in populating the levels, the ortho- 
to-para ratio will be significantly smaller than 3 (Martini 
et al. 11997)1 . Hence the ortho-to-para ratio can be used as 
a diagnostic to determine whether FUV fluorescence plays 
a role in exciting the ro-vibrational states. 

To calculate the ortho-to-para ratio from our observed 
line fluxes, we used Eq. 5 from Hoban et al. (1991) with 
the transition probabilities and level energies taken from 
Martin et al. 1)1996)1 . Since we observed many more lines 
than they did, we chose to calculate the ortho-to-para ra- 
tio from each combination of ortho and para states that 
was observed in our spectrum. We also assigned an uncer- 
tainty to each ratio, based on the estimated uncertainties 
in the line fluxes and calibration. We then calculated the 
weighted average of all ratios, resulting in an ortho-to-para 
ratio of 2.77 ± 0.19. Note that the uncertainty quoted here 
is the uncertainty in the average value, not the spread of 
the individual measurements. This value is in good agree- 
ment with the expected ratio of 3 for collisionally excited 
molecular hydrogen. 
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5.3. The (1,0)S(1)/(2,1)S(1) ratio 

A diagnostic that has been commonly used to discrim- 
inate between shocks and FUV-pumped fluorescence, is 
the (1,0)S(1)/(2,1)S(1) ratio. In shocks this ratio is typ- 
ically ~10 (e.g., Shull & Hollenbach 119781 and references 
therein) while for pure radiative fluorescence it is 1.8 
(Black & Dalearno I1976[l . This ratio is essentially based 
on the same physics as the comparison between the ro- 
tational and vibrational temperatures, but is a useful di- 
agnostic when no other H2 lines have been observed. The 
dereddened value we obtain for this ratio is 16.2, thus in- 
dicating shock excitation. 

5.4. Thermal excitation 

Based on the evidence presented in the previous subsec- 
tions, we can state that the H2 emission is collisionally 
excited by other atoms or molecules in the gas. These col- 
liding particles could obtain their energy either from ab- 
sorbing stellar UV photons with energies less than 11.2 eV 
(more energetic photons cannot be present, they would 
dissociate H2), or from shocks. In order to investigate this 
point further, we examined the 12 000 K post-AGB star 
model from run 4 with graphite dust and no dust forma- 
tion in the post-AGB wind (van Hoof et al. I1997[) . This 
model should give a fair approximation for the physical 
conditions in IRAS 16594—4656. Using the distance es- 
timate from Sect. and the known angular dimensions 
of IRAS 16594—4656, we noticed that the radial sizes in 
this model were approximately a factor of two too large. 
Hence we reduced all radii by a factor of two and increased 
all densities by a factor of four. We then recalculated the 
model using Cloudy 96 beta 4 (Ferland[2002 ) . This showed 
that in regions where molecular hydrogen is abundant, the 
electron temperature is 210 K or lower. This value is sig- 
nificantly lower than the rotational temperature of 1440 K 
which we derived above. It is therefore very unlikely that 
H2 is thermally excited by UV heated gas. 

6. Discussion 

Based on the evidence presented in the previous section, 
we can conclude that the H2 emission is collisionally ex- 
cited by shocks in IRAS 16594-4656. 

Molecular shocks generally can be described as "C- 
shocks" or as "J-shocks" , depending on their shock veloc- 
ity v s , the ambient ionization fraction and the strength 
of the component of the ambient magnetic field perpen- 
dicular to v a (Draine 1980). A J-shock is a shock where 
the heat deposition length is short compared to the cool- 
ing length. J-shocks are sufficiently powerful to dissociate 
molecules. C-shocks occur at relatively low velocities and 
moderate to low ionizations and depend upon the pres- 
ence of a magnetic field. In C-shocks, gas is accelerated 
and heated by collisions between charged particles and 
neutral particles. C-shocks cool mainly through molecular 
material and grains. In C-shocks the radiation is emitted 



as the gas is being heated in the shock front; in J-shocks it 
is emitted after the impulsive heating event, downstream 
behind the shock front. 

The H2 emission in IRAS 16594—4656 appears strong: 
much stronger than Br7 for instance. Molecular lines pro- 
duced by collisional excitation of H2 will be strongest 
if the collisions are not energetic enough to dissoci- 
ate H2 and lower its abundance (Hollenbach & McKee 
1989). Such conditions exist in C-shocks because the gas 
is heated gradually and remains molecular. Hence large 
columns of material with temperatures between 1000 K 
and 3000 K can be produced which generate copious 
amounts of H2 emission. Therefore the strong H2 emis- 
sion in IRAS 16594—4656 argues in favor of a C-shock. 

We compared our line ratios of the H2 emission with 
the intensities of ro-vibrational transitions of H2 from C- 
shock models by Le Bourlot et al. ( 2002). Our values are 
consistent with a 20 - 30 kms -1 C-shock propagating in 
10 3 cm -3 material. They are not consistent with any of 
their other models where the C-shocks have either a higher 
velocity or impinge on higher density gas. 

However, C-type shocks depend upon the presence of 
a magnetic field. There is no direct evidence for the pres- 
ence of a magnetic field in IRAS 16594-4656. Su et al. 
(2003J observed 10% polarization close to the central star, 
which makes it likely that the EDE harbors one. Direct 
detection of magnetic fields is very difficult, and has only 
been attempted a handful of times. Nevertheless, a detec- 
tion has been achieved in four post-AGB stars and PNe. 
Miranda et al. I|2001(l detected a toroidal magnetic field 
in the EDE of the young PN K3-35. Greaves lj2T)TI!2)l de- 
tected a toroidal magnetic field around NGC 7027 with 
an effective scale of 5000 AU and correctly oriented to 
collimate the observed bipolar winds, and of a dimension 
similar to the base of the outflow. They also detected a 
magnetic field in CRL 2688. Bains et al. il2l)l)3i measured 
for the first time the magnetic field strength in a post- 
AGB star. They found B = +4.6 mG at 1612 MHz and 
B = +2.5 mG at 1667 MH in OH 17.7-2.0. 

The lack of significant ionization (no radio emission or 
lines typical of H 11 regions have been detected) is consis- 
tent with the prerequisites of a C-type shock. 

However, metastable lines such as [0 1] 6300 A and 
[C 1] 9850 A have been detected. These are only produced 
in J-shocks. C-shocks, which have large molecular abun- 
dances, tend to rapidly convert any pre-existing atomic 
oxygen into molecular form (Hollenbach & McKee 1989). 
The observation of ionic lines such as [Fen] are also a 
strong indication of J-shocks, if one has eliminated the 
possibility of an H 11 region as the origin. Few ions exist in 
C-shocks, which only form when the shocked gas is largely 
neutral and molecular. 

Another indication in favor of the presence of a J-shock 
is the wind speed of 126 kms -1 derived from the Ha P- 
Cygni profile (Van de Steene et al. l2000bjl . This velocity is 
too high for a C-shock, but consistent with velocities for 
a J- type shock (Hollenbach & McKee [THEH • 
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[Ci] 9850 A, [Ni] 1.040 /jm, and [Fell] 1.257 & 
1.644 /im are all predicted to be strong near- infrared 
emission lines in J-type shocks, while in the optical [N ii] 
6548 k 6584 A, and [S n] 6716 & 6731 A would be strong. 
Comparing the detected near-infrared line fluxes to the 
model predictions presented in Fig. 8 of Hollenbach & 
McKee l|1989jl . we find that they are compatible with an 
80 kms -1 J-shock impinging upon 10 3 cm -3 material, or 
alternatively a 70 kms -1 shock impinging upon 10 4 cm' 3 
material. This would be consistent with a 126 kins -1 wind 
colliding against the side walls of an hourglass nebula at 
an angle of approximately 45° w.r.t. the normal, assuming 
that the side walls are expanding at a rate of 16 kms -1 
(Loup et al. 1990). This value for the impact angle is in 
fair agreement with the wide opening angle of the lobes 
seen in the image in Su et al. |2001). 

In the optical (Reyniers, private communication) and 
infrared spectra no forbidden lines from ionized species 
other than [Fe n] were detected. A weak detection of [N n] 
6548 & 6584 A would be expected in the optical spectrum 
according to the 10 3 cm -3 model of Hollenbach & McKee 
(1989). For a given shock velocity, the ionized states will 
become relatively weaker for higher densities since the 
overall degree of ionization will shift downwards. Taken 
at face value this could point to a relatively high density 
(at least 10 4 cm -3 ) at the edge of the lobes. However, 
according to the 10 4 cm -3 model, a weak detection of 
[Ni] 1.040 /im would be expected, which was not observed 
in our infrared spectrum either. The critical density of 
this line at T c = 10 4 K is approximately 10 7 cm" 3 , using 
data from Zeippen l|1982[> and Dopita et al. I|1976[l . This 
value is larger than for any of the other lines listed above, 
ruling out an extremely high density and collisional de- 
excitation of the line as the cause of the non-detection. 
Consequently the non-detections of [Ni] 1.040 /im and 
[Nil] 6548 & 6584 A seem to contradict each other and 
could be indicative that IRAS 16594—4656 has a low ni- 
trogen abundance. Bipolar PNe with high C/O and low 
N/O ratios are known (e.g., J 900, Perinotto 1991), al- 
though bipolar PNe are usually associated with type I PN, 
which have high nitrogen abundances. One should keep in 
mind though that the Hollenbach & McKee Ijl989|l mod- 
els use solar abundances. The carbon-rich chemistry in 
IRAS 16594—4656 could have an important effect on the 
models since it will alter the cooling significantly. Deeper 
observations are needed to settle this point. 

^From the previous discussion it is seems that slower 
C-type and faster J-type shocks both are present in 
IRAS 16594—4656. Because the ionization potential of 
neutral iron is 7.87 eV, and the dissociation energy of H2 
is 4.48 eV (Graham et al. 1987), in principle H2 and Fe + 
cannot coexist in the same region in substantial quanti- 
ties. Hence, the C- and J-shocks must occur in different 
regions in the nebula. 

C-shocks would occur where the shock velocity is low, 
hence where the fast wind impinges almost tangentially 
upon the circumstcllar material. The H2 emission would 
originate from this region. We postulate that the C-shocks 



occur where the fast wind is funncled through the EDE 
of IRAS 16594—4656. This is in agreement with Guerrero 
et al. l|2000l) who noted that, for PNe, in all cases where 
the emission observed in the H2 line is more than 5 times 
greater than the Br7 emission, the H2 emission is mainly 
concentrated in a bright ring, the "waist" of the bipolar 
nebula. For IRAS 16594—4656 the Br7 line is much weaker 
than the H 2 (1,0) S(l) line. The ratio H 2 (1,0) S(l) to Br 7 
is 8.4 after correction for extinction. As mentioned above, 
the EDE is likely to harbor a magnetic field, which is a 
prerequisite for a C-shock. 

However, molecules have been known to survive in 
dense condensations and cometary knots (e.g., Speck et 
al.|2D02). Hence it cannot be excluded that some H2 emis- 
sion is present in the lobes, excited by shocks caused by a 
molecular outflow impinging on the AGB envelope (e.g., 
the Egg Nebula, Cox et al. [2000). The large difference 
seen between the IRSPEC and SOFI value of H 2 (1,0) S(l) 
may support this. However, no collimated outflow nor high 
density blobs have been observed in IRAS 16594-4656 yet. 

J-shocks occur where the shock velocity is high. A first 
possibility is that the [Ci] and [Fell] line emission origi- 
nates in a J-shock along the walls of the bipolar lobes, 
where the fast wind impinges with a high transverse ve- 
locity. This is in agreement with the observations of Hora 
& Latter lfT59"fi| and Welch et al. fiMty who imaged 
the shock excited [Fell] line emission in the young PN 
Hubble 12. The shape of this PN resembles a tilted hour- 
glass. The [Fen] emission is located along the side walls 
of the hourglass in this object and its morphology differs 
strikingly from that of the H2 emission, which is mostly 
dominated by UV excitation. The geometry of this PN is 
very similar to IRAS 16594-4656, but Hubble 12 is a PN 
and hence at a much more evolved stadium. 
A second possibility is that the [Fe 11] emission originates 
much closer to the central star, in the post-AGB wind 
itself. Shock waves induced by stellar pulsations were pro- 
posed to explain the [Fen] emission in Mira variables 
(Richter et al.[2003). Post-AGB stars are of course hotter 
than Miras, but they are usually variable and still pulsat- 
ing. Fokin et al. (200Q argue that stellar pulsations are 
forming shocks in the atmosphere of the post-AGB star 
HD 56126. Miras are believed to have very high mass loss 
rates, but the mass loss rate in post-AGB stars may still be 
quite large (up to 10™ 5 -10" 6 M Q yr" 1 ,Gauba et al.|HS|), 
as may be indicated by the strong P-Cygni Balmer lines 
in the optical spectrum of IRAS 16594—4656. Moreover, 
post-AGB stars have higher wind velocities, which lead to 
the same normalized [Fe 11] peak fluxes at lower pre-shock 
densities. 

A third possibility is that the J-shocks would occur where 
matter transferred from a binary companion hits an accre- 
tion disk. At this stage there is no corroborating evidence 
that either a binary companion or an accretion disk exists 
in IRAS 16594-4656. 

High resolution images in H2 and in selected 
metastable lines will clarify which shock mechanism occurs 
where, while high spectral and spatial resolution infrared 
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spectra will clarify the kinematics around the central star 
and in the bipolar lobes. 

7. Conclusions 

We examined the near-infrared spectrum of 
IRAS 16594—4656. It shows strong H2 emission lines and 
some typical metastable shock excited lines such as [Fe 11] 
1.257 & 1.644 jim. We argue that the molecular hydrogen 
emission is mainly collisionally excited. Its strength indi- 
cates that the H2 emission originates in C-type shocks. 
However, the metastable lines, and especially the [Fell] 
emission lines, indicate the presence of J-type shocks. 
These shocks don't usually coexist. Hence we postulate 
that the H2 emission originates mainly where the stellar 
wind is funneled through an EDE and therefore impinges 
almost tangentially upon the circumstellar material. 
The [Fen] emission either occurs along the walls of the 
bipolar lobes where the transverse shock velocity would 
be higher, or could originate much closer to the central 
star in shocks in the post-AGB wind itself, or possibly in 
an accretion disk. 

IRAS 16594-4656 is unique in the sense that it allows 
us to study shock excited lines which arc not compromised 
by ionization or fluorescence. Further high resolution near- 
infrared spectra are currently being obtained in order to 
clarify the kinematics of this object and the shock mech- 
anisms. 
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